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Abstract 

A global model of a slowly rotating perfect fluid ball in general relativ- 
ity is presented. To second order in the rotation parameter, the junction 
surface is an ellipsoidal cylinder. The interior is given by a limiting case 
of the Wahlquist solution, and the vacuum region is not asymptotically 
flat. The impossibility of joining an asymptotically flat vacuum region 
has been shown in a preceding work. 

1 Introduction 

In a preceding paper |Q , henceforth called I, the impossibility of matching the 
Wahlquist metric to an asymptotically flat vacuum domain was shown. This 
result is not too surprising in the light of investigations in (5) , where the tendency 
of the matching conditions to be overdetermined has been pointed out. However, 
it would be very embarrasing from the point of view of general relativity if this 
matching turned out to be impossible to any vacuum region. In this paper, the 
problem of matching of the slowly rotating Wahlquist metric to a more general 
vacuum exterior is investigated to a precision of quadratic order in the angular 
velocity. 

An approximation scheme for slow rotation was first introduced by Brill and 
Cohen || || to first order in the small angular velocity. In our paper we are 
using the the formalism developed by HartleQ, taking into account quadratic- 
order terms in the power series expansion in the angular velocity f2 of the fluid. 
The metric of both regions has the form 
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ds 2 = (1 + 2h)A 2 dt 2 - (1 + 2m)B 2 dr 2 

-(1 + 2k)C 2 \dd 2 + sin 2 i? [dp + (fi - w) rfi] 2 } (1) 

where the static, spherically symmetric state is described by the functions A, B 
and C depending only on the radial coordinate r, while the functions to, h, m 
and k can, in general, depend both on r and The rotation potential to is of 
first order in the angular velocity f2, and the functions h, m and k are of order 

n 2 . 

In the Wahlquist interior domain, the rotation potential to is a function of 
the radial coordinate alone. Hence, from the juction conditions it follows, that 
to is independent of $ in the vacuum region as well. In our model of the space- 
time we drop the condition of asymptotic flatness, and we perform the matching 
with the most general vacuum metric with to — to{r). Likewise the expansions 
of the second-order metric functions in Legendre polynomials are sought in the 
form h = ho + /12-P2 (cos m = mo + m?.Pi (cost?), k — fc2P2(cosi9), where the 
functions ho, /12, fiOj m 2 and k 2 depend only on the radial coordinate r, and 
P 2 (costf) = (3cos 2 i9- l)/2. 

In Sec. 2 of this paper, we present the form of these functions for the 
vacuum domain, and investigate the effect of the non asymptotically flat part 
of the perturbations on the curvature. In Sec. 3., the junction conditions are 
calculated, and the constants determining the exterior metric as functions of 
the unperturbed radius T\ are given in full detail. 



2 The Vacuum Exterior 

In this section we consider the form of the vacuum metric to the required ac- 
curacy. The unperturbed metric is described by the Schwarzschild solution, 
A 2 = 1/B 2 = 1 — 2M/r and C = r. The solution for the perturbed metric is of 
the form (Q), with O = and 

2aM 

u = -78- • P) 

where the additive constant of integration has been removed by a rigid rotation. 
Integration of the second-order metric functions yields 

1 fa 2 M 2 r \ 
1 f a2M2 \ 



2M 



2M \ „ M 



h 2 = 3cir (2M - r) log I 1 j+a 2 — (M + r) 

M , 9 9N r-M ( 2M \ , 

+2ci — 3r 2 - 6Mr - 2M 2 ) — +1 r 2 Ql 5 

r 2M — r \ r J 
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2M\ „M . 



3ci(r 2 - 2M 2 ) log ( 1 - ^ ] - a 2 ^-(2Af + r) 
M 



2ci — (2M 2 - 3Mr - 3r 2 ) + (2Af 2 - r 2 ) gi (6) 



m 2 = ba — hi . (7) 



v 



In this approximation, the slowly rotating solution is characterized by the mass 
M, the first order small rotation parameter a, and the second order small con- 
stants ci, C2 and q\. 

With the special choice gi = 0, the metric is asymptotically flat, and this 
form was used in [[[J to show the impossibility of matching the Wahlquist solution 
to an asymptotically flat vacuum. Since the asymptotic behaviour changes 
completely when qi becomes nonzero, the far field behaviour cannot be treated 
in a perturbative way, and the present series expansion is certain to hold only 
within an open neighborhood of the junction surface. This is also indicated 
by the fact that the terms containing the integration constant q\ tend to r 2 as 
r — > oo. 

To explore the effects of the qi perturbations on the curvature, we use a 
canonical locally nonrotating Lorentz tetrad for the metric (□): 



r 2 „ . o A 1 d 



e = {l-h+^stf*)^ (8) 
ei - (l-m)A^- (9) 

62 = (1 - fc)r_1 ^ (10) 

In this tetrad, the gravitoelectric (Ei) and gravitomagnetic (Hi) fields are 
defined in terms of the Riemann tensor as follows H : 

Ei = i?oioi E% = R0202 E3 — i?oio2 (12) 

Hi = i?0123 H2 = — i?0213 H3 — i?0223 • (13) 

In the asymptotic region, r — > oo, the gravimagnetic part of the curvature goes 
to zero. The gravielectric components tend to the finite values 



lim Ei = (l — 3 cos 2 i?) gi (14) 

lim E 2 = (l-3sin 2 tf)gi (15) 

'— *CX3 

lim E% — 3 sin cos $ qi . (16) 



The algebraic structure of the Weyl tensor at infinity, (^4j)-(|l6|), guarantees that 
the components cannot all be transformed to zero, and hence the space-time 
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cannot be asymptotically flat. However, these values are obtained by means 
of perturbative calculations, and higher-order corrections may contribute by 
divergent terms]] to the limiting values (p^)-(16). 



The angular behavior of the gravielectric field indicates that a quadrupolar 
mass distribution at large distance may act as the source of the deviations from 
asymptotic flatness. To show this, let us consider a large sphere with radius R, 
and with a non-uniform surface density distribution 



X>2^(cost?) , (17) 



i=0 



where \xn are constants. In the weak-field approximation, linearizing around 
the background Minkowski metric and expanding in powers of 1 /R, the metric 
near the center of the sphere can be written as 



where 



ds 2 = (1 + 2iP)dt 2 - (1 - 2iP){dr 2 + r 2 d$ 2 + r 2 sin 2 ■ddtp 2 ) (18) 



^ = -\^R-^2r 2 P2{cos§)+o(jp) (:i!>) 



(With the conventions used in I the gravitational constant is g^.) The terms 
can be absorbed by rescaling the t and r coordinates, while the /12 terms have 
exactly the same angular and radial dependence as the qi terms in the far field 
exterior vacuum region in (||)-((7|). The higher than quadrupole mass distribu- 
tion coefficients can be arbitrary, which shows that they cannot be determined 
from a slow rotation formalism which takes into account only up to quadratic 
order terms in the angular velocity Q. Instead of using the linearized gravity 
approximation, an exact static but not spherically symmetric Weyl-class vac- 
uum metric could also be constructed. But since for slow enough rotation there 
is always a region where the r 2 q\ terms are small whereas the other terms in 
the metric are negligible, no new insight would result from that analysis. 

In a general stationary axisymmetric vacuum exterior solution the suitable 
hypersurfaces for matching to an interior fluid region are determined in the 
paper of Roos[Q. In the limit of no rotation, the matching surface is the history 
of the sphere r = T\. For slow rotation the surface becomes an ellipsoid, and its 
deformation is described by 

r = n +a 2 [xo + X2P2 (cos 1?)] (20) 

where a is the small rotational parameter and xo arL d X2 are constants to be 
determined by the matching conditions. The expressions for the normal vector 
and the extrinsic curvature of the matching surface have the same forms as in 
I. 



1 We remind the reader that the asymptotic region may lie outside the domain of conver- 
gence of the power series expansion in the angular velocity. 



4 



3 Matching with the Wahlquist solution 



In I we have computed the second order form of the Wahlquist metric H in 
Hartle's coordinates. The metric functions and the zero-pressure matching sur- 
face is expressed in I in terms of the constants (Iq and n characterizing the static 
configuration and in terms of the the small parameter tq which is proportional 
to the angular velocity. These results will be used now for matching with the 
vacuum solution of the previous section. The radial coordinate in the fluid re- 
gion is denoted by x instead of r. In the limit of no rotation the fluid region is 
described by the Whittaker metric H, and the matching surface is the sphere 
characterized by x = x%, where x\ cotxi = n 2 (cf. I). 

For slow rotation the matching surface S is an ellipsoidal cylinder character- 
ized by a vanishing pressure and by the embedding condition (po|). We equate 
with each other the respective induced extrinsic curvatures K(y) and Krw) of 5, 
in the vacuum and in the Wahlquist region. Hence, the equations of matching 
are 

dsf v) \ s = ds 2 (w) \ s K (v) \s = K {W) \ S . (21) 

where ds 2 \s is the induced metric. The values of the metric coefficients and 
their derivatives on S are given by a power series expansion in tq in the fluid 
and in a in the vacuum regions, respectively. We apply a rigid rotation in the 
fluid region by setting ip — > ip + fit where f2 is a constant. Then we re-scale the 
interior time coordinate t by 

t ^ c 4 {l + r 2 c 3 )t (22) 

with further constants C3 and C4 to be determined from the matching conditions. 

Substituting in the matching conditions (|2l|), and Taylor expanding to sec- 
ond order in the angular velocity, we get a set of linear equations for the pa- 
rameters xoj X2, ci, C2, C3 and qi . Here we list the solution of the matching 
equations for the constant C3 of the interior time scaling and for the constants 
Qi, ci, C2, Xo and Xi °f the vacuum domain in terms of the radius x\ of the 
Whittaker fluid ball, the rotation parameter a given by 

rn 2x\ cos 2 x\ — 3 sinxi cosxi + x\ . , 

a = (23) 

3 cos x\ sin x\ cos x\ — x\ 

and the density /iq. The results of I for the matching to zero order in the angular 
velocity 

n = ^ QXiro (24) 

Gsinxi cosxi 

(25) 



hold without any change, 



M = — -^(k 2 -cos 2 xi ) (26) 

2 l/2 

ri = -T^-sinx! (27) 

c 4 = cosxi . (28) 
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The second-order matching conditions have the solution 

COSXlito r/ , _ 2 „ • 2 \ 3 

c 3 = -T-TTo 5 : (15cos Xi + 7sin Xi)xJ 

72Zg COS" 1 Xi SIMi 

— (18cos 2 Xi + 4sin 2 x\)x\ sinxi cosxi + 12cos 3 x\ sin 3 xi (29) 

— 3xi cos 2 xi sin 4 xi — 9xi cos 4 xi sin 2 xi] 



x 



C2 = 



9 9 9 3 

a cos Xi/igxf 



2^/3 s i n 3 Xl 

[51x 3 cos 3 xi — 27x 2 cos 4 x\ sinxi — 27xi cos 5 xi sin 2 xi 



+ 3 cos 4 xi sin 3 x\ — 24xi cos 3 X\ sin 4 xi 



34x 2 cos 2 xi sin 3 xi + 8x 3 cos x\ sin 2 xi — x\ sin xi] 
3a 2 li cos 5 xi/XgX 5 
Z2 Z3 sin 3 xi 



+ 2 4 - 3 log(smxiCosxi/xi) (30) 



a 2 li cos 5 xi/x 2 x 5 
Z 2 /! 4 sin 3 xi 

27 1/9 • 1/2 1/2 
— a 13 cos ' xi sin ' xi/i 



*2 *3 sin xi 



25/2/22-3/2 

x [9x 3 cos 2 xi - 27xi cos 4 x\ sin 2 xi (32) 
+ 18 cos 3 x\ sin 3 xi — 9xi cos 2 x\ sin 4 x\ + 2x\ cos xi sin 3 x\ + x 3 sin 2 xi] 

3? 3 cos 3 / 2 xi^o /2 2 4 

Yn = tt\ 7 fx, cosxi — 2xi cos X\ sin x\ (33) 

23/2xl /2 ; 2 2 sin 3 / 2 xi 

+ cos 3 xi sin 2 xi + X! sin 3 xi — cosx! sin 4 x x ) 

— 3/3 COS 3 / 2 X\llrJ 2 2 . ,33 „ 3 .2 

y 2 = T7K 7-^ (x-i smxi + 4xV cos xi — 8x1 cos xi sm xi (34) 

23/2x} /2 ; 2 sin 3 / 2 x 1 x 1 111 

+ 3 cos 2 X! sin 3 x x — 4x! cosxi sin 4 Xi) 
where 

11 = 3x 2 + 6x1 cosxi sinxi — cos 2 xi sin 2 xi — 8 sin 2 xi (35) 

12 = xi(3cos 2 xi + sin 2 xi) — 3sinxi cosxi (36) 
h = xi — sinxi cosxi (37) 

It is easy to show that qi is non-zero in the allowed range of Xi, < Xi < |. 

4 Conclusions 

In this paper we have shown that there exist configurations where a Wahlquist 
fluid ball in an ambient vacuum domain is kept in equilibrium. In our model 
space-time the components of the curvature tensor tend to finite values at in- 
finity. A far-away quadrupole mass distribution can ensure the required non- 
asymptotically flat nature of the exterior vacuum. The solution is linearization 
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stable in the sense that there exists a family of exact solutions corresponding 
to our approximate solution, which is valid in the fluid and in at least an open 
vacuum region surrounding it. This follows from the work of Roos H, [l^ show- 
ing the existence and uniqueness in a neighbourhood of the matching surface in 
the general axisymmetric and stationary case. 
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